Fas/FasL of pacific cod mediated apoptosis.

Fas and Fas ligand (FasL) pathway plays important roles in virus defense and cell apoptosis. In our previous work, nervous necrosis virus (NNV) was discovered in Pacific cod (Gadus macrocephalus), and the Fas ligand (PcFasL) was up-regulated when NNV outbreak, however, signal transmission of Fas/FasL in fish are still unclear. In the present study, Pacific cod Fas (PcFas), PcFasL and Fas-associating protein with a novel death domain (PcFADD) were characterized. The predicted protein of PcFas, PcFasL and PcFADD includes 333 aa, 90 aa and 93 aa, separately. 3-D models of PcFas, PcFasL and PcFADD were well constructed based on reported templates, respectively, even though the sequence homology with other fish is very low. The transcript levels of PcFas increased gradually from 15 day-post hatching (dph) to 75dph. PcFas was significantly up-regulated when cod larvae had NNV symptoms at 24dph, 37dph, 46dph, 69dph, and 77dph. Subcellular localization revealed that PcFasL was located in the cytoplasm, while PcFas was mainly located in the cell membrane. Exogenous expressed PcFasL of 900 µg/mL could kill the Epithelioma papulosum cyprinid (EPC) cells by MTT test, but low concentration has no effect on the cells. qPCR analysis showed that overexpression of PcFas could significantly up-regulate the expression of genes related to Fas/FasL signaling pathway, including bcl-2, bax, and RIP3, while overexpression of PcFasL significantly up-regulate the expression of caspase-3, caspase-9, and MLKL. Overexpression of PcFas or PcFasL could induce EPC apoptosis significantly by flow cytometry, which was consistent with the results of caspase-3 mRNA level increasing. The results indicated that NNV could induce apoptosis through Fas/FasL signal pathway.

Expression and antibacterial analysis of galectin-8 and -9 genes in mandarin fish, Siniperca chuatsi.

Galectin-8 and galectin-9 belong to tandem repeat-type galectins, and in the present study, these two genes were cloned in mandarin fish Siniperca chuatsi. The open reading frame (ORF) of the mandarin fish galectin-8 and galectin-9 contains 942, and 1008 bp, encoding 313 and 335 amino acids, respectively. As a conserved feature, an N-terminal carbohydrate recognition domain (CRD), and a C-terminal CRD were observed in each of the two galectins in mandarin fish. In healthy fish, galectin-8 and -9 were constitutively expressed in all organs/tissues examined, and their expression can be induced following the stimulation of LPS and poly(I:C). It is obvious that galectin-8 had a higher increase at mRNA level following the stimulation of poly(I:C). It is further demonstrated that mandarin fish galectin-8 inhibited the growth of Flavobacterium columnare and Streptococcus agalactiae, and in addition to the two species of bacteria, galectin-9 inhibited also the growth of Edwardsiella piscicida, which provides the basis for further understanding their antibacterial role in immune response of mandarin fish.

Identification of differentially expressed genes in gills of tiger puffer (Takifugu rubripes) in response to low-salinity stress.

Salinity is an important abiotic factor for aquatic organisms. In fish, changes in salinity affect physiological responses and alter the immune system. Takifugu rubripes is an important economic marine fish, and mechanisms of T. rubripes adaptation to salinity changes need to be further documented. In this study, a transcriptome sequencing technique was used to analyse genes that were differentially expressed in the T. rubripes gill after low-salinity stress for 30 d, and differential gene expression was further validated by quantitative real-time PCR (qPCR). After assembly, 385 differentially expressed genes (DEGs) were identified, including 182 upregulated genes and 203 downregulated genes. The DEGs were assigned to Gene Ontology (GO) classes with a total of 1647 functional terms. Most DEGs were assigned to biological process (984; 59.8%) followed by molecular function (445; 27.0%) and cellular component (218; 13.2%). Further KEGG analysis allocated 385 DEGs to 95 KEGG pathways. After q-value correction, 7 pathways (Glycolysis/Gluconeogenesis; Biosynthesis of amino acids; Carbon metabolism; Fructose and mannose metabolism; Pentose phosphate pathway; Metabolism of xenobiotics by cytochrome P450; and Glycine, serine and threonine metabolism) remained significant. qPCR results indicated that the transcripts of six selected genes sharply increased after 30 d of low-salinity stress. Low-salinity stress obviously increased SLC39A6, SLC5A9, NKAα1, CYP1A1, CYP1B1, and GSTA expression. In contrast, the genes encoding Aldoaa, GPI, FBP2 and GAPDH exhibited downregulation. In addition, three solute carrier (SLC) genes selected from the DEGs were further studied for differential expression patterns after low-salinity exposure, and the results showed that the SLCs were upregulated in T. rubripes after 72 h of low-salinity exposure. This investigation provides data for understanding the molecular mechanisms of fish responses to low-salinity stress and provides a reference for rationally setting salinity levels in aquaculture.

Conserved function of Pacific cod Caspase-3 in apoptosis.

Apoptosis plays a significant role in the cellular immune responses against infections, especially those related to viruses. Across various species, Caspase 3 is a prominent mediator of apoptosis and participates in the cell death signaling cascade. However, its role remains relatively unknown in cod fish. In this study, we aimed to reveal the role of Pacific cod Caspase-3 (GmCasp3) in apoptosis and its evolutionary position. Our results showed that the GmCasp3 cDNA contains an open reading frame of 864 nucleotides; that codes for 287 amino acids long protein with a molecular weight of 32.03 kDa. The sequence alignments and 3-D model indicated that GmCasp3 contained highly conserved domains, such as "QACRG", "GSWFI" and "HG" active sites, however, the phylogenetic tree analysis revealed that both GmCasp3 and Atlantic cod caspase-3 clustered together are far from the high vertebrate branch, indicating they are at a lower position in vertebrate evolution. Red fluorescent labeling vector pDsRed2-C1-GmCasp3 was constructed and it was transfected into EPC cell lines. The result showed that GmCasp3 protein was distribute in the protoplasm and expressed in apoptotic cell debris. Moreover, the GmCasp3 enzyme activity increased with the increased post-transfection analysis time, while the genome DNA was visibly fragmented at 36 h post transfection. Flow cytometry analysis showed that the proportion of apoptosis cells increased from 12 h to 24 h post transfection. In conclusion, the conserved functions of GmCasp3 in apoptosis indicated that Pacific cod has the similar apoptotic characteristics as other animals.

Transcriptome from Pacific cod liver reveals types of apolipoproteins and expression analysis of AFP-IV, structural analogue with mammalian ApoA-I.

Apolipoproteins (Apos), transporting the lipids through the lymphatic and circulatory systems, are associated with kinds of diseases. Additionally, type IV antifreeze protein (AFP-IV) was related evolutionarily with apolipoproteins. However, the information of Apos in fish was limited. In this study, ApoA-I, ApoA-I-2, ApoA-IV, Apo E, ApoB-100-like and AFP-IV were sequenced from Pacific cod (Gadus macrocephalus) liver transcriptome using Illumina HiSeq 2000, and their 3-D models were constructed based on the most confidence templates ever reported in mammals. Interestingly, the model of G. macrocephalus AFP-IV, named GmAFPIV, is quite similar to the structure of ApoA-I. GmAFPIV includes 689 bases with a complete open reading frame encoding 125 amino acids. Sequence alignment of GmAFPIV showed 30% to 50% similarity with that of other species except Gadus sp. Expression levels of GmAFPIV were found in a decreasing manner in liver, intestine, gill, brain and gonad. Heterologously expression of the GmAFPIV protein was expressed in Escherichia coli and purified to immunize New Zealand rabbits. The survivors of E. coli in 60 µg/mL of GmAFPIV are more than that in the 30 µg/mL group after stored in -20 °C and -80 °C, indicating high concentration of GmAFPIV could protect E. coli avoiding the damage from ice crystal. The subcellular localization of GmAFPIV showed that the green fluorescence was mainly observed in the cytoplasm, indicating GmAFPIV play roles in the cytoplasm. It was concluded that GmAFPIV may function not only as an antifreeze protein but also as an apolipoprotein transporting lipids in fish.

Molecular characterization of caspase members and expression response to Nervous Necrosis Virus outbreak in Pacific cod.

Multiple functions of caspases include normal cell turnover, proper development and function of the immune system, hormone-dependent atrophy, embryonic development and chemical-induced cell injury. During artificial propagation of Pacific cod, Gadus macrocephalus, high mortality occurred during early development stages. Here, we performed various analyses on the cDNA and protein sequences of six different G. macrocephalus caspases namely GmCasp3, 6, 7, 8, 9 and 10, and tried to investigate the contributions of caspase family to the development and Nervous Necrosis Virus (NNV) resistance. Sequence analysis of GmCaspase proteins showed that each caspase shared conserved domains like "HG", "QACXG (X for R, G or Q)" and "GSWF" except GmCasp10. Sequence alignment and phylogenetic tree showed that GmCasp8 and GmCasp10 were quite different from those of other fishes. 3-D models indicated that structure of GmCasp3 is very conservative, but GmCasp6, 7, 8, 9 and 10 are less conservative. Tissue distribution analysis showed that six Gmcaspases mRNA transcripts were detected in tissues of intestine, gill, thymus, head-kidney and spleen with different abundance, but Gmcasp7 were not detected in the brain. GmCasp3 transcript was kept at very low level in the early development stages, while the expression levels of GmCasp6, 7, 8, 10 were different at various development stages. GmCasp8 level seemed to be much higher than other caspases in the heads of 65dph and 75dph juveniles. To understand the role of caspases during NNV outbreak, modulation in expression of each Gmcaspases were investigated. The results showed that GmCasp3 transcript level increased significantly when NNV broke out, while GmCasp7, 8, 9 and 10 in cod heads decreased obviously at 69dph and 77dph. The results suggest that caspases in Pacific cod should be diverse in their structure and function, and their unique features and response to NNV outbreak add more evidences for the specificity of immune system in Pacific cod.

Digital gene expression analysis of Takifugu rubripes brain after acute hypoxia exposure using next-generation sequencing.

The adverse effects of hypoxia are confined to biochemical, physiological, developmental and behavioral processes, especially injury of the brain. In this study, a subset of genes in the brain of Takifugu rubripes were analyzed using digital gene expression (DGE) profiles and next-generation sequencing after acute hypoxia. Among 32 differentially expressed genes, 29 were up-regulated and 3 were down-regulated following hypoxia exposure. Using Gene Ontology analysis, it was found that transcription and translation, metabolism, and the stress response were affected by exposure to hypoxia. KEGG analysis revealed that the neuroactive ligand-receptor interaction pathway was significantly enriched in hypoxia-exposed T. rubripes. To further confirm the differential expression of genes, quantitative real-time PCR was performed to test six candidate genes, with the following five genes exhibiting the same expression patterns as the sequencing results: Proto-oncogene c-fos, Kruppel-like factor 2, immediate early response 2, proopiomelanocortin A and rhodopsin. This work is the first to identify and annotate genes in T. rubripes affected by hypoxia stress. This investigation provides data for understanding the molecular mechanism of fish adaptation to hypoxia and provides a reference for rationally setting dissolved oxygen levels in aquaculture.

Evidence for and characterization of nervous necrosis virus infection in Pacific cod (Gadus macrocephalus).

A mortality rate higher than 90% was observed in a larva-rearing facility for Pacific cod, Gadus macrocephalus, in China. Larvae showing clinical signs of infection were collected. Initial suspicion of nervous necrosis virus (NNV) infection was confirmed by sequencing, absolute quantification real-time PCR (A-qPCR), and electron microscopy. The nucleotide sequence of RNA2 was 1,375 bases long (GenBank no. KM576685), coding for a single ORF corresponding to the capsid protein from residues 21 to 1034. Phylogenetic analysis of the capsid protein sequence showed that PCNNV belongs to the barfin flounder NNV (BFNVV) genotype. An amino acid sequence alignment revealed 39 differences between the cold- and warm-resistant viral groups, suggesting that PCNNV evolved under temperature selection. The 3-D structure of the predicted capsid protein was modeled to identify potential epitopes, and the gene was expressed in Escherichia coli, yielding a protein with a molecular mass of 55 kDa. During PCNNV outbreaks, the viral copy number was found to reach 10(7) per ng of total RNA, which could be considered the lethal copy number of NNV in cod. The gonads, eggs, fertilized eggs and asymptomatic cod fry were all positive for PCNNV, indicating viral vertical transmission as the main source of the viral load. The amount of virus in the apparent healthy fry or survivors seemed to decrease gradually with development. These results might lead to efficient diagnostic methods to help farmers select NNV-free broodfish for cod breeding.

Transcriptomic analysis and biomarkers (Rag1 and Igµ) for probing the immune system development in Pacific cod, Gadus macrocephalus.

Mortality (>90%) is a big concern in larval rearing facilities of Pacific cod, Gadus macrocephalus, limiting its culture presently still in the experimental stages. Understanding the immune system development of G. macrocephalus is crucial to optimize the aquaculture of this species, to improve the use of economic resources and to avoid abuse of antibiotics. For the transcriptome analysis, using an Illumina sequencing platform, 61,775,698 raw reads were acquired. After a de novo assembly, 77,561 unigenes were obtained. We have classified functionally these transcripts by Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG). 27 genes mainly related to hematopoietic or lymphoid organ development and somatic diversification of immune receptors have been reported for the first time in Pacific cod, and 14 Ig heavy chain (µ chain) locuses were assembled using Trinity. Based on our previous achievement, we have chosen Rag1 and Igµ as immune system development biomarkers. Full length cDNA of Rag1 and Igµ as biomarkers were obtained respectively using RACE PCR. Concerning Rag1, the deduced amino acid of Rag1 and protein immunodetection revealed a Rag1 isoform of 69 kDa, significantly different from other fish orthologs, such as Oncorhynchus mykiss (121 kDa). Phylogenetic analysis reveals a unique immune system for the Gadus genre, not exclusive for Atlantic cod, among vertebrates. Meanwhile, full length cDNA of Igµ included an ORF of 1710 bp and the deduced amino acid was composed of a leader peptide, a variable domain, CH1, CH2, Hinge, CH3, CH4 and C-terminus, which was in accordance with most teleost. Absolute quantification PCR revealed that significant expression of Rag1 appeared earlier than Igµ, 61 and 95 dph compared to 95 dph, respectively. Here we report the first transcriptomic analysis of G. macrocephalus as the starting point for genetic research on immune system development towards improving the Pacific cod aquaculture.

Differential gene expression profile of the calanoid copepod, Pseudodiaptomus annandalei, in response to nickel exposure.

To better understand the underlying mechanisms of reactions of copepods exposed to elevated level of nickel, the suppression subtractive hybridization (SSH) was used to elucidate the response of the copepod Pseudodiaptomus annandalei to nickel exposure at the gene level. P. annandale is one of a few copepod species that can be cultured relatively easy under laboratory condition, and it is considered to be a potential model species for toxicity study. In the present study, P. annandalei were exposed to nickel at a concentration of 8.86 mgL(-1) for 24h, after which the RNA was prepared for SSH using unexposed P. annandalei as drivers. A total of 474 clones on the middle scale in the SSH library were sequenced. Among these genes, 129 potential functional genes were recognized based on the BLAST searches in NCBI and Uniprot databases. These genes were then categorized into nine groups in association with different biological processes using AmiGO against the Gene Ontology database. Of the 129 genes, 127 translatable DNA sequences were predicted to be proteins, and the putative amino acid sequences were searched for conserved domains (CD) and proteins using the CD-Search service and BLASTp. Among 129 genes, 119 (92.2%) were annotated to be involved in different biological processes, while 10 genes (7.8%) were classified as an unknown-function gene group. To further confirm the up-regulation of differentially expressed genes, the quantitative real time PCR were performed to test eight randomly selected genes, in which five of them, i.e. α-tubulin, ribosomal protein L13, ferritin, separase and Myohemerythrin-1, exhibited clear up-regulation after nickel exposure. In addition, MnSOD was further studied for the differential expression pattern after nickel exposure and the results showed that MnSOD had a time- and dose-dependent expression pattern in the copepod after nickel exposure. To the best of our knowledge, this is the first attempt to investigate the toxicity effects of nickel on a copepod at molecular level.

Characterization of RAG1 and IgM (mu chain) marking development of the immune system in red-spotted grouper (Epinephelus akaara).

In vertebrates, lymphoid-specific recombinase protein encoded by recombination-activating genes (RAG1/2) plays a key role in V(D)J recombination of the T-cell receptor and B-cell receptor. In this study, both RAG1 and the immunoglobulin M (IgM) mu chain were cloned to characterize their potential role in the immune defense at developmental stages of red-spotted grouper, Epinephelus akaara. The open reading frame (ORF) of E. akaara RAG1 included 2778 nucleotide residues encoding a putative protein of 925 amino acids, while the ORF of the IgM mu chain had 1734 nucleotide residues encoding 578 amino acids including variable (VH) and constant (CH1-CH2-CH3-CH4) regions. E. akaara RAG1 was composed of a zinc-binding dimerization domain (ZDD) with a RING finger and zinc finger A (ZFA) in the non-core region and a nonamer-binding region (NBR), with a zinc finger B (ZFB), the central and C-terminal domains in the core region. Tridimensional models of the ZDD and NBR of E. akaara RAG1 were constructed for the first time in fishes, while a 3D model of the E. akaara IgM mu chain was also clarified. The RAG1 mRNA was only detected in the thymus and kidney of 4-month and 1.5-year old groupers using qPCR, and the RAG1 protein was confirmed using western blotting and immunohistochemistry. The IgM mu mRNA was examined in most tissues except the gonad. RAG1 and IgM mu gene expression were observed at 15 dph (days post-hatching) and 23 dph respectively, and increased to a higher level at 37 dph. In addition, this was the first time that the morphology of the E. akaara thymus was characterized. The oval-shaped thymus of 4-month old fish was clearly seen and there were amounts of T lymphocytes present. The results suggested that the immune action of E. akaara was likely to start to develop around 15 dph to 29 dph. The transcript level of the RAG1 gene and the number of lymphocytes in the thymus between 4-month and 1.5-year old groupers indicated that age-related thymic atrophy also occurs in fishes. The similar functional structures of RAG1 and IgM protein between fish and mammals indicated that teleost species share a similar mechanism of V(D)J recombination with higher vertebrates.